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The imminent detection of gravitational waves will trigger precision tests of gravity through 
observations of quasinormal ringing of black holes. While General Relativity predicts just two po- 
larizations of gravitational waves, the so-called plus and cross polarizations, numerous alternative 
theories of gravity predict up to six different polarizations which will potentially be observed in 
current and future generations of gravitational wave detectors. Bekenstein's Tensor- Vector-Scalar 
(TeVeS) theory and its generalization fall into one such class of theory that predict the full gamut 
of six polarizations of gravitational waves. In this paper we begin the study of quasinormal modes 
(QNMs) in TeVeS by studying perturbations of the scalar field in a spherically symmetric back- 
ground. We show that, at least in the case where superluminal propagation of perturbations is 
not present, black holes are generically stable to this kind of perturbation. We also make a unique 
prediction that, as the limit of the various coupling parameters of the theory tend to zero, the QNM 
spectrum tends to 1/V2 times the QNM spectrum induced by scalar perturbations of a Schwarzschild 
black hole in General Relativity due to the intrinsic presence of the background vector field. We 
further show that the QNM spectrum does not vary significantly from this value for small values of 
the theory's coupling parameters, however can vary by as much as a few percent for larger, but still 
physically relevant parameters. 

PACS numbers: 04.50.Kd, 04.30.-w, 04.70. Bw, 04.80. Cc 



I. INTRODUCTION 



Perturbations of black holes have been intensely stud- 
ied during the past decades with relation to black hole 
stability, astrophysical implications of gravitational wave 
detection and, more recently, gauge/gravity dualities (see 
the reviews PHI])- From an astrophysical perspective, 
the impending observation of gravitational waves from 
black holes will herald new and unprecedented tests of 
gravity through two main avenues. The first is through 
the possible detection of gravitational wave polarizations 
that are not predicted by General Relativity (GR) such 
as scalar "breathing modes" and extra tensor and vector 
degrees of freedom 1 . The second method is a direct probe 
of the no-hair theorem of GR through observations of the 
quasinormal ringing of black holes known as quasinormal 
modes (QNMs). General Relativity predicts that the si- 
multaneous measurement of the frequency and damping 
time of a single QNM is sufficient to determine both the 
mass and angular momentum of a black hole [p34l6| . 
implying the detection of more than one QNM from a 
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1 See H forpossiblc gravitational modes in metric theories of 
gravity, and 0,@| as well as [9l 4llH and references therein for dis- 
cussions of possible detections with inteferometer and spherical 
gravitational wave detectors resepctivcly. 



single black hole can give rigorous constraints on the no- 
hair theorem - one of the major science goals for ground 
based and space based gravitational wave detectors [e.g. 
ITfjl - fli^ . An alternative approach to the measurement of 
the frequencies of black hole QNMs based on a connection 
between the emission of gravitational waves and strong 
gravitational lensing has recently been elucidated in Ste- 
fanov et al. [2fjj . 

A recently popular alternative theory of gravity that 
predicts both extra polarizations of gravitational wave 
propagation as well as a violation of the no-hair theorem 
is Bekenstein's Tensor- Vector- Scalar (TeVeS) theory [2l[ 
and its generalization [22|-[24|. In this paper we take a 
first step towards exploring the perturbations of black 
holes within TeVeS and its generalization by studying 
perturbations of the background scalar field. This en- 
ables us to look at both the stability of black holes in the 
theories, as well as looking at how the QNM frequencies 
and damping times arc affected by the various parame- 
ters of the theories. 

TeVeS is a covariant generalization of the Modified 
Newtonian Dynamics (MoND) paradigm [25|, which at- 
tempts to explain the discrepancy between observed and 
predicted mass distributions on galactic scales. As a tool 
for interpreting galactic scale observations, MoND has 
proven extremely successful (see [26j for a review) . How- 
ever, as a non-covariant theory it is ill-equipped to ex- 
plain the full gamut of astrophysical and cosmological 
observations. To this end, TeVeS has recently been em- 
ployed in an attempt to explain, amongst other things, 
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gravitational lensing [27H3il . the cosmic microwave back- 
ground power spectrum j32l l33j , large-scale cluster sur- 
veys [H[ and even type la supernova [H, [3(| , while re- 
taining the successful predictions of MoND in the weak 
acceleration limit (see |37| for a recent review of TeVeS). 

In the Newtonian regime, TeVeS has been shown to re- 
produce the parametrized post-Newtonian coefficients to 
a level consistent with solar system experiments [2l| . In 
the strong field regime however, comparatively little has 
been studied. Giannios 38] first solved the field equa- 
tions for static, spherically symmetric, vacuum space- 
times within TeVeS. He found two branches of solutions 
dependent on the degrees of freedom allowed in the vec- 
tor field. Giannios' solution was, however, plagued by the 
fact that the scalar field was necessarily negative for some 
radii, implying superluminal propagation of scalar waves 
[2lj . Sagi and Bekenstein [3^] remedied this situation 
by showing that another branch of solutions to the field 
equations exist that has the same physical metric, but 
where the scalar field is positive throughout the space- 
time. Sagi and Bekenstein [39[ also found the charged 
Reissner-Nordstrom solution within TeVeS and provided 
a detailed study of black hole thermodynamics. Lasky 
ct al. [!(| then found the Tolman-Oppcnheimer-Volkoff 
solution in order to study neutron star structure, work 
which has since been extended to include slow rotation 
jilf and fluid & spacetime perturbations [H, EH • 

Meanwhile, trouble was brewing with the ori gina l for- 
mulation of the TeVeS field equations. Scifert 4J] first 
showed that the Schwarzschild-TeVeS solution is unsta- 
ble to linear perturbations for experimentally and phc- 
nomenologically valid values of the various coupling pa- 
rameters. Contaldi et al. [22j then showed that the vec- 
tor field is prone to the formation of caustics in a variety 
of simple dynamical situations, analogously to Einstcin- 
T/thcr theories where the vector field is described by a 
Maxwellian action. Finally, Sagi [lij showed that the 
vector field is constrained by the cosmological value of 
the scalar field in such a way that it prevents the scalar 
field from evolving. 

To overcome the aforementioned issues, Contaldi et al. 
H2, Skordis [23| and Sagi 24| independently provided a 
generalization of TeVeS by taking a more general action 
for the vector field which is motivated by Einstein- jEthcr 
theory. Complicating an already complicated theory is 
a complicated process. It has implied that little rigor- 
ous work has been achieved looking at the structure of 
the field equations. Skordis [23| studied the cosmological 
equations of the theory and showed they are identical to 
the equations governing the original TeVeS cosmology up 
to a rescaling of Hubble's constant. Sagi [H| looked at 
the parametrized post-Newtonian (PPN) parameters of 
the theory and found them not to be in conflict with so- 
lar system experiments for suitable values of the coupling 
parameters. 

It should be mentioned that the generalized TeVeS 
theory is also not without its problems. In two papers, 
Mavromatos et al. 45 1 and Ferreras et al. [46| have shown 



that observations of gravitational lensing and galactic ro- 
tation curves provide an inconsistent parameter space for 
the theory. Although their results were based on the 
original TeVeS theory, their formulation of the problem 
carries over to the generalized theory based on their use 
of a vector field with only a temporal component (the 
proof of this is the same as that given in Lasky (47| in 
the strong-field regime). However, it is exactly this point 
that diminishes the robustness of their claims. One main 
driving force for introducing the vector field into the the- 
ory is to provide sufficient gravitational lensing without 
requiring dark matter (2lj . It is therefore not surpris- 
ing that suppressing the degrees of freedom of this vec- 
tor field by setting all spatial components to zero also 
suppresses the degree to which gravitational lensing ob- 
servations can be made to be consistent with other ob- 
servations. With the extra complexity of the vector field 
in the generalized theory, combined with the relatively 
primitive state of parameter space estimation, one would 
expect that gravitational lensing observations could be 
induced into conforming with galactic rotation curve ob- 
servations with a more robust form of the vector field. 

The study of gravitational wave propagation in Gen- 
eralized TeVeS was first specifically broached by Sagi 
[is]] , who studied linear perturbations of the field equa- 
tions to look at the speed and form of gravitational wave 
propagation. Sagi showed that there exist six different 
modes of gravitational waves (as opposed to two in GR) 
which propagate at four different speeds, all different to 
the speed of light. As expected, the various propaga- 
tion speeds of the gravitational wave modes were found 
to be dependent on the various coupling parameters in 
the theory. It is interesting to note that this is seem- 
ingly in contrast with the result of Kahya and Woodard 
[49j and Desai et al. who predict gravitational waves 
will propagate at the same speed independently of the 
coupling parameters. The obvious discrepancy between 
these two results implies further investigation is neces- 
sary. 

Spherically symmetric vacuum, charged and perfect 
fluid solutions of the gen eralized TeVeS field equations 
were studied in Lasky [47| where the vector field was as- 
sumed to contain only a temporal component. Under 
these symmetry assumptions it was shown that these so- 
lutions are the same as for the original TeVeS theory by 
Giannios [38j], Sagi and Bekenstein [39[ and Lasky et al. 
[ioj respectively, where only a rescaling of the vector field 
coupling parameters is required. That is, given the afore- 
mentioned rescaling, the form of the background scalar, 
vector and tensor fields, and therefore also the form of 
the physical metric, are exactly the same in the two the- 
ories. In Lasky [47| it was further shown that this result 
is not gencralizable - solutions with time dependence or 
more complicated geometry in the vector or tensor fields 
will be different between the two theories. As we are 
dealing in the present paper with perturbations of the 
spherically symmetric vacuum solution, we are dealing 
with the solution of both the original and generalized 
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versions of TeVeS. Moreover, as the difference between 
these two theories is in the vector field, the perturbations 
we perform of the scalar field also hold for both theories. 
That is, all work presented in this paper is applicable 
to both the original version and generalized version of 
TeVeS. As such, unless explicitly mentioned, we herein 
refer to TeVeS to mean the all-encompassing generalized 
version of the theory. 

TeVeS admits the Schwarzschild solution as a possible 
geometry [3^, H3] , and most likely the Kerr solution as is 
the case with a majority of alternative theories of grav- 
ity [5l[ (although we note there exists various theories of 
gravity that do not admit the Kerr solution, for example 
Chern-Simons gravity |H3|). Therefore, electromagnetic 
observations of the spacetime surrounding black holes 
may not allow for the distinction between GR, TeVeS and 
alternative theories. For example, observing the dynam- 
ics of stars orbiting close to Sgr A* will soon be yielding 
fruitful information about the spacetime structure of that 
black hole [Ejl,|54j|. This problem however, is well approx- 
imated by the stars acting like zero-mass test particles, 
implying one is only probing geodesies of the background 
spacetime. Therefore, measuring the spacetime to be a 
Kerr black hole using this method does not rule out alter- 
native theories of gravity that admit Kerr as a possible 
geometry. The perturbations of such spacetimes (rele- 
vant for gravitational waves) however, depends on the 
specific field equations of the theory, and hence will dif- 
fer between theories. This implies QNMs of various Kerr 
(or Schwarzschild) geometries behave differently and can 
be used to provide rigorous tests of the theory of gravity. 

The article is set out as follows; in section [TTI we pro- 
vide a brief primer on the relevant field equations for the 
generalized TeVeS theory, reviewing the structure of the 
background spacetime in section IIII1 In section IIVI we 
derive the wave equation governing scalar field perturba- 
tions and also discuss analytic results in the limit of small 
couplings. In section IV A[ we analyse the stability of the 
black hole solutions and then in section IV Bl wc compute 
the QNM spectrum. Throughout the article Greek in- 
dices range from ... 3 and antisymmemetrization and 
symmctrization are respectively denoted by square and 



The modified Einstein field equations are 
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II. GENERALIZED TEVES EQUATIONS 

TeVeS is built upon three dynamical fields; the Ein- 
stein metric, g^, a time- like normalized vector field, A 11 , 
and a dynamical scalar field, (p. A physical metric, g^, 
in which clocks and rulers are measured, is related to the 
other three fields according to 



^(g^+A^-e^A^. (1) 



This implies that variations in the tensor, vector or scalar 
fields can be measured through their coupling to the 
physical metric. 
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where G^ is the Einstein tensor associated with the Ein- 
stein frame, is the physical stress-energy tensor and 
t m „ and 0^„ are the effective stress-energy terms asso- 
ciated with the scalar and vector fields respectively. In 
particular 
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where k is the scalar field coupling constant, /i, a function 
of the theory's coupling parameters and the scalar field, 
is associated with the MoND acceleration scale, T (fi) is 
related to the interpolation function in MoND and hence 
is not a priori predicted by the theory and £ is a fixed 
length scale associated with the free function. For the re- 
mainder of this article we only consider the strong-field 
regime of the theory, whereby [i = 1 is an excellent ap- 
proximation (for more details see [U [M, H^]). In 
this case, Contaldi et al. [22j showed that the function T 
logarithmically diverges, but is exactly cancelled in the 
field equations, implying this function has zero contribu- 
tion in this limit. Moreover, it has been shown (39| that 
this limit of the strong-field approximation, whereby we 
allow fj, = 1, is correct out to at least a million times the 
gravitational radius off a black hole. Not only is this well 
into the asymptotically flat region of the black hole, but 
for the purpose of QNM calculations this is more than 
sufficient for the external boundary condition. That is, 
it is sufficient in an astrophysical sense, that when we are 
treating QNMs and their boundary conditions at "infin- 
ity" , we are still dealing with the strong-field limit of the 
theory. 

One key difference between the original version of 
TeVeS developed by Bekenste in J21| and the current ver- 
sion with Maxwellian action 22, 231 is the introduction 



of three extra vector coupling constants which describe 
the relative strengths of the individual terms in the vec- 
tor field action. The four vector field coupling constants 
are denoted by K, K +1 K2 and K4, where the origi- 
nal theory is regained by setting the last of these three 
parameters to zero. The vector field contribution to the 
effective stress-energy in the modified Einstein equations, 
MI/ , can now be expressed as a sum of terms associated 
with each of the vector field coupling constants, plus a 
term associated with the Lagrange multiplier, A (which 
ensures normalization of the vector field); 
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where 
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in TeVeS. Giannios [38[ first solved the original TeVeS 
equations of Bckenstein for the case of a spherically sym- 
metric, static, vacuum spacetime. He showed that there 
exist two branches of solutions based on the form of the 
vector field; one branch where the vector field is tempo- 
rally aligned and another when the vector field has both 
a temporal and radial component. Lasky (47j then found 
the equivalent solution in the case of a purely temporal 
vector field for the more general TeVeS theory. In that 
work it was shown that the solution in the generalized 
TeVeS theory is the same as that in the original TeVeS 
theory with a simple substitution of the parameters of 
the theory 

Following the notation of Giannios [38|, we write the 
physical line element in isotropic coordinates in terms of 
two constants, r c and a; 
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Here, F M „ := W[ V A^, S MV := V^^) and A^ := 
A a S/ a A fl . 

Variation of the total action with respect to the vector 
field gives the vector field equation 

K\7 a FK a + K+V a S afl + K 2 W (V a A a ) + XA" 
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Contraction of this equation with the vector field isolates 
the Lagrange multiplier. This subsequent equation can 
then be used in the modified Einstein equation, in par- 
ticular in the term expressed in equation ([§]), such that 
the system is fully determined. 

The final field equation is that of the scalar field 
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It is of interest for the current work how the TeVeS field 
equations limit to the GR equations. This is achieved by 
continuously limiting k and all of the K's to zero, as 
well as taking t — > oo. This was shown in detail for the 
original TeVeS theory by Bekenstein [2l[ , his sections III 
C and D. 



III. SPHERICALLY SYMMETRIC, STATIC, 
VACUUM SOLUTIONS 



In order to discuss perturbations of black holes, it is 
worth spending some time discussing the current sta- 
tus of spherically symmetric, static, vacuum solutions 



where dfl 2 :— d8 2 + sin 2 9 dip 2 . The two constants in 
the physical metric are related to the various coupling 
parameters of the theory, including the "scalar mass", 
, as well as the characteristic gravitational radius, r g , 
according to 



(13) 
(14) 
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where JC := K + K + — K4. Note that all quantities in 
equation (|14[) are greater than or equal to zero, implying 
a > 0. Moreover, there are only two situations where a 
can be identically zero. Firstly, if both the gravitational 
radius and the scalar field coupling constant, r g and fc, 
both vanish. Not only does this reduce the theory to the 
Einstein- iEther theory, but it also implies that the metric 
becomes Minkowski. The second possibility is that r c 
diverges, which is unphysical. It is therefore only possible 
that a is strictly greater than zero. 

Equation (fT2"|) describes more than just a black hole 
spacetime. Giannios [38[ showed that r c is a black hole 
event horizon if and only if a = 2. This is based on 
the requirement that the surface at r = r c must have a 
finite surface area, and also the singularity must be re- 
movable. Evaluating the physical Ricci scalar shows that 
it diverges for values of a < 2 and also 2 < a < 4. Mean- 
while, as the surface area is proportional to the g rr (r c ) 
component of the metric, one finds that this diverges for 
a > 2. One therefore has a black hole solution only when 
a = 2. When < a < 2 and 2 < a < 4 one has a 
naked singularity at r = r c , and a > 4, r c represents a 
removable singularity with a divergent surface area. 

As mentioned, the only solution given by the metric 
(fl"2"|) that represents a black hole is that with a = 2. One 
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immediately notes that this is exactly equivalent to the 
Schwarzschild metric in isotropic coordinates. 



tion IV Bl that perturbations in TeVeS have a different 
spectrum. 



A. Equivalence with Brans-Dicke Theory 



B. Vector and Scalar Field Parameters 



It is further interesting to note that the metric of equa- 
tion (|12[) is exactly equivalent to the Brans 'type F metric 
[55| in Brans-Dicke scalar tensor theory [56| (with the as- 
sociation a = 2Q and x = from the notation of Scheel 
et al. [H?}) [HI]. Note that in the Brans I solution, setting 
X = necessarily implies that the parameter Q is unity, 
and hence the Schwarzschild solution is recovered. This 
is because the Brans I solution has an extra algebraic 
equation which links Q and x- 

A popular representation of the Brans type I metric 
can be found by taking the coordinate transformation 
given by 



R{r) = r(l + r c / r y 



(15) 



from the metric (|T2j) and defining M := 2r c , one gets 
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In particular, this is the line element expressed as equa- 
tion (7) of Campanelli and Lousto [59| with the associa- 
tion m = —n = a/2 — 1. For a majority of the article we 
will deal with the case of a = 2, for which the line element 
(fT6| becomes exactly the Schwarzschild line element. 

The association with the Brans-Dicke theory is inter- 
esting given the vast complexity of the field equations in 
TeVeS as compared with Brans-Dicke theory. However, 
the association is due to the simplifying assumption that 
the vector field is purely timclike. It was shown for the 
original version of TeVeS that black hole solutions with 
a non-zero radial component for the vector field produce 
significantly more complicated geometries than the case 
where the radial component vanishes [38j . This complex- 
ity will only increase when one analyses the black hole 
solutions with non-zero radial vector fields in the fully 
general theory. 

It is clear that, while the background spacetime in 
TeVeS and Brans-Dicke theories are equivalent, pertur- 
bations of these spacetimcs will differ as the scalar field 
equation of TeVeS is not a Klein-Gordon equation but 
includes contributions from the background vector field. 
In fact, whilst scalar field perturbations of Schwarzschild- 
Brans-Dicke black holes [equation (fl6)) with a = 2] have 
the same QNM spectrum as that of scalar perturbations 
of Schwarzshild black holes in GR 2 we shall show in scc- 



The above metric (|12j) is that found by Giannios |38(. 
He displayed trepidation towards this solution because 
the scalar field becomes negative close to the horizon, 
which therefore allows for supcrlumin al p ropagation of 
scalar waves [2l[ . Sagi and Bekenstein [39| overcame this 
issue by showing that Giannios [38| overlooked a branch 
of solutions where the scalar field was everywhere posi- 
tive, although they did this only for the case where a = 2. 
Using equations (|T5|) and (|T?|) we can re-derive the extra 
solution of Sagi and Bekenstein [39| with the generaliza- 
tion that a remains arbitrary. In this case, one can show 
that the scalar field throughout the spacetime is given by 



ip(r) = if c + S± In 
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where <p c is the cosmological value of the scalar field, 
which can be determined through cosmological observa- 
tions, and 5± can be found to be 
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Note here that when a — 2, equation (|T8f reduces to 
exactly equation (67) of Ref. [39[ (where K. = K as 
they were working in the original TeVeS which has K + = 
K 2 = K 4 = 0). 

Given that < K, < 2 and k > 0, one can trivially 
show that (5+ > 0. Moreover, after some extra work a 
conditions is derived for 5- that is 



a 2 < 



2-K 



S_ < 



(19) 



which further implies that 5- < for all a < 2. A 
majority of the results of this paper are for the cases 
where a = 2, which is the only black hole solution of 
the field equations. In this case it is clear that <5_ is 
always negative, a fact that will be important for the 
perturbation analysis in section fVl 

It is interesting to note that there exists a special case 
of the above functions whereby <5_ vanishes identically, 
implying that ip(r) = tp c throughout the spacetime. It is 



2 This was not explicitly shown in Kwon et al. |6C| who analyzed 



the stability of such perturbations, however the wave equation 
they derive is equivalent to the Regge- Wheeler wave equation for 
scalar perturbations of Schwarzschild black holes in GR. 
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apparent that this is only the case where a 2 = 8/(2 — /C) 
(or when /C = which is not of interest for the present 
work). Whilst this case exhibits a constant scalar field 
throughout the spacetime, there still exists an irremov- 
able naked singularity at r = r c . Although the stability 
of such objects is of great interest, we shall see below 
that the wave equation governing the perturbations for 
this case is no simpler than the generic case when there 
is a naked singularity at the horizon. As such, we leave 
this for further exploration in a latter paper. 

The generic behaviour of the scalar field throughout 
the spacetime is of great importance to the present dis- 
cussion due to the presence of superluminal propaga- 
tion of scalar waves. Indeed, Bekenstein [2l[ showed 
that, in the eikonal approximation, scalar perturbations 
in the Newtonian limit of the theory travel with veloc- 
ity v = exp (— 2ip) /v2- This implies that scalar per- 
turbations travel superluminally if and only if ip < 0. 
Equation (fTTj) shows that the scalar field necessarily di- 
verges as r — > r c . When S± < 0, the scalar field diverges 
to positive infinity, and indeed the scalar field is every- 
where positive (given ip c > 0). However, when 5± > 0, 
the scalar field diverges to negative infinity as one ap- 
proaches the critical radius, r c , implying for some values 
of the radial coordinate the scalar field is negative. 

Finally, for the remaining sections we will require the 
form of the vector field in the prescribed background 
spacetime, which can be shown to be 
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We note that the various coupling parameters of the the- 
ory are mainly expressed here in the e v term [sec equation 
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IV. PERTURBATION EQUATIONS 



We express the perturbed scalar field as 

Sip 

ip = ip + Sip, where — << 1, (21) 

where ip and Sip are respectively the background and 
perturbed scalar fields. Since the scalar field pertur- 
bations decouple from perturbations of the vector and 
tensor fields [48|, we are free to let perturbations of the 
vector and tensor fields vanish, i.e. SA^ = Sg^ v = 0, 
throughout the article, and we leave the analysis of these 
perturbations to a future article. 

The perturbation of equation ([!]) implies 



-2Sp 
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Given this equation, together with the fact that the phys- 
ical metric of the background spacetime is diagonal, this 



can be expressed component-wise as 



Sg tt = 28ipg t t 



and 
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From the above form of the perturbed physical met- 
ric it is obvious that a perturbation of the scalar field 
only affects the diagonal components of the line ele- 
ment. These are the so-called breathing modes of gravita- 
tional waves which are not present in GR. The detection 
of these modes is possible using current intcrfcrometric 
techniques, although requires multiple suitably oriented 
detectors to be operating simultaneously (see Ref. [1] 
and references therein), spherical resonant-mass detec- 
tors (for example (^4llL loll l62j). or through the differ- 
ence in expected energies between the standard plus and 
cross polarizations in standard intcrfcrometric detection. 
The detection of such modes would instantly falsify GR, 
however not detecting these modes will simply act to con- 
strain both sources and parameter spaces in various the- 
ories. 

The scalar field equation ([lip is simple to perturb, and 
turns out to be 



[(g a0 - A a A?) d a Sp] = 0. 



(24) 



For the remainder of the article we work only with the 
black hole case in which the metric parameter a = 2. In 
the Appendix we show the more general form of the wave 
equation with arbitrary a, however we leave the study of 
the stability of such spacetimes to future work. 



A. The wave equation 

Performing the standard seperation of variables in 
Schwarzschild coordinates 



S(p (t,R,t 



*(t,R) 
R 



EE 

1=0 m=—i 



Y, 



ira i 



(25) 



where the Yg m are the standard spherical harmonics, 
leads to the wave equation 



-2e 4 ^ 



<9 2 $ <9 2 $ 



dt 2 OR 2 



-11- 



2M 
~R~ 



1(1 + 1) 
R 2 



2M 



$ = 0, 
(26) 



where the standard tortoise coordinate, R+(R), has been 
defined according to 



dR 2M 
dRZ ~ ~~ ~R~' 



(27) 



An important aspect of this wave equation (|26[) is that it 
is almost identical to the wave equation governing scalar 
perturbations of a Schwarzschild black hole in GR, the 
difference being the factor of 2 exp (4ip) preceding the 
second-order time derivative. In the weak field limit (i.e. 
taking large R) equation (|26[) reduces to a wave equation 
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where the velocity of the perturbations is exp (— 2ip) / \/2, 
in agreement with the velocity of scalar perturbations de- 
rived by Bekenstein [2l| and Sagi |48j . 

Consider now the limit of equation (|26p as the back- 
ground scalar field goes to zero. This limit is not well 
defined in the sense that the background scalar field nec- 
essarily diverges as i? — > 2M for 5± ^ 0. However, if we 
take the limit as K,, k and <p c tend towards zero, we find 
the background scalar field tends to zero for all values of 
R 2M. In this limit we see that the coefficient of the 
d 2 /dt 2 term tends to two, which remains distinct from 
the general relativistic limit. Consider now a rescaling of 
the temporal coordinate t = \f2t' , which acts to reduce 
the wave equation to exactly the Regge- Wheeler equa- 
tion for scalar perturbations of Schwarzschild black holes 
in GR. If we now follow standard methods and assume a 
harmonic time dependence for the perturbations 



JLU) s t' 



(28) 



where u> and lu s are the QNMs for our TeVeS black hole 
and the GR Schwarzschild black hole respectively, then 
we find the following simple relation 



lim 



V2 



(29) 



That is, as the theory's various coupling parameters tend 
to zero, the QNM frequencies and damping times are re- 
lated to the Schwarzschild QNM frequencies and damp- 
ing times through a factor of \[2. We note again that 
the GR Schwarzschild black hole spectrum is equivalent 
to the QNM spectrum of Brans-Dicke black holes [60j |. 
and hence the TeVeS QNMs also deviate by the same 
relation to the Brans-Dicke QNMs. We have confirmed 
the result of equation ([2Ut numerically in section IV B I 

Naively, one may expect QNMs in TeVeS to tend ex- 
actly to the QNMs of GR in the limit as the theory's 



parameters tend to zero. However, this expectation 
breaks down as the perturbation equation being consid- 
ered herein does not exist in GR because of the lack of a 
vector field, and hence the derivation of the wave equa- 
tion is not applicable. 



B. Wave equation in Isotropic Coordinates 

Equation (|26|) is not appropriate for studying details 
of QNMs due to the non-trivial coefficient preceding the 
time derivative. Interestingly, performing a coordinate 
transformation from the Schwarzschild coordinates to the 
original isotropic coordinates defined by the transforma- 
tion (fl"5)). together with the substitution 



*(i,r) 



M 
1+ 2^) 



(30) 



leads to the standard form of the wave equation 
d 2 d 2 

V ] * = 0. 



dt 2 dr 2 

Here, we have introduced a new coordinate defined by 



(31) 



dr \ 2r 



3-2S± 



1 - 



M 
27 



-l+25± 



(32) 



which is a "tortoise-like" coordinate, in the sense that it 
maps r € (M/2, oo) to G (—00,00), providing —1 + 
28± < 0. Note that in isotropic coordinates the horizon 
occurs at r = M/2. 

The potential, y(r+), is now given by the complicated 
expression 



V 



3 -4<p c 

2r 2 



-2(3-25±) 

27 J 



M v 2(l-2«±) r 



M 
27 



2(l-fa)-^ 
r 



M 
1+ 27 



M 



(33) 



Throughout section |V] we will look for QNMs of the 
above system of equations. For these purposes it is con- 
venient to assume a harmonic time dependence 

V(t,r) = §(r)c^*, (34) 

where uj = wr c + iwi m , which implies equation pip be- 



comes 

+ (. 2 -V)* = O. (35) 

In the above definition, wr c represents the frequency of 
the QNM while u!\ m is its damping time. Any mode with 
w im < will grow exponentially in time, and hence rep- 
resents an unstable mode. 
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C. Behaviour of the potential 

Determining stability and finding the QNM spectrum 
is now a one-dimensional scattering problem that is de- 
pendent solely on the form of the potential, V, defined in 
equation (|33|). It is therefore worth spending some time 
establishing various properties of this function. 

We firstly note that the potential vanishes at spacelike 
infinity, independent of the value of ip c . The behaviour 
approaching the horizon, r = M/2, is more difficult as 
this is dependent on the value of 5±. Bearing in mind 
that, for the black hole case, 6- < and 6+ > 0, we find 



0.03- 



0.02- 

V 

0.01 - 




1 23456789 10 11 



r/2M 



lim V(r) 

. M + 



for 5- 
— oo for <5+ 



(36) 



As discussed, the case of S + necessarily has regions of 
the spacetime that allows for superluminal propagation 
of scalar field perturbations. While it has been heavily 
debated whether this is physically allowed to occur in na- 
ture (see for example Ellis et al. [63|, Bruneton (64| and 
references therein), it is clear that this induces a pathol- 
ogy on the wave equation governing scalar perturbations 
in the form of a divergent potential function. When we 
are discussing such superluminal perturbations we must 
be careful with our definition of "horizon" . Indeed in 
this case the surface at r = M/2 is a photon horizon 
associated with the null cones of the physical metric. Su- 
perluminal propagation of scalar field perturbations im- 
plies that a perturbation at some r < M/2 could escape 
through the horizon at r = M/2 out to infinity. We 
believe that this induces the pathological behaviour of 
the potential for the S+ case, however we reserve further 
exploration of the effect of this to a future article. 

It is also of interest to determine when the potential 
function is negative as such regions of a potential can rep- 
resent bound states which imply growing modes. Con- 
sider, for the moment, the £ = case. From equation 
(f3"3")l . one can show that for 6- the potential is negative 
in the region 



M M 

— < r < , 

2 " 2 (!-<*_)' 



and for 5+ the potential is negative for 



M Mil- <L 

—— < r < ^— — - 

2 2 



(37) 



(38) 



It is trivial to show that for all values of £ there exists 
a negative portion of the potential, however the largest 
negative region is realized for the £ = case described 
above. A generic plot of the potential is shown in figure 
[TJ for the S- case with £ = 0, 1 and 2. We note that at 
this scale the difference in behaviour between the S- and 
S + cases can not be seen. However, when one zooms into 
the region close to the horizon, figure [3J the difference in 
the respective functions is apparent. 



FIG. 1: Potential, V, for the scalar perturbation wave equa- 
tion for a black hole (i.e. a = 2) with K, — 0.1, k = 0.01 
and ip c = 0.003. Note that r is an isotropic radial coordinate 
implying the horizon of the spacetime is at r = M/2. This is 
the case of <5_ , although at this range plots of 5+ look similar 
(the effect of the sign choice in S± becomes more relevant for 
higher ^'s). Here, the thick black, dashed blue and dotted red 
lines represent £ = 0, 1 and 2 respectively. 




1 .002 1 .004 

r/2M 



1.006 



1.008 



FIG. 2: A zoomed in view of figure [T] for 5- (top panel) 
and S+ (bottom panel). As stated in the text, as r -> M/2 
the potential goes to zero for all £'s for the 5- case, however 
diverges to negative infinity for the case of 8+. This divergent 
behaviour is associated with the superluminal propagation of 
the scalar perturbations. 



V. RESULTS 
A. Stability analysis 

The presence of a single unstable mode can do severe 
damage to a theory as a viable alternative theory of grav- 
ity. The stability analysis of spherically symmetric black 
holes is generally rendered rather simple due to the prop- 
erty that any growing, unstable mode, must necessarily 
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be purely imaginary (Konoplya and Zhidcnko [for exam- 
ple see |4j). That is, if an unstable mode exists (i.e. such 
that u>i m < 0), it must have zero real part, and therefore 
will not be oscillatory in nature. 

We have performed time evolutions of equation ([26]) 
for various values of the spherical harmonic I. The I = 
mode has the largest negative region in the potential, and 
hence is the most likely to be unstable. In figure [3J we 
plot the time evolution of a typical I = mode (the figure 
uses the example where JC = k = 0.1 and ip c = 0.003). 
Consistently with the case with scalar perturbations of 
the Schwarzschild black hole in GR, the £ = mode has a 
large imaginary part (see section l"VB[) . implying the sig- 
nal is significantly damped after only a few oscillations. 
Whilst this makes it difficult to extract robust measure- 
ments of the frequencies and damping times for these 
modes, one can clearly see from figure [3J that oscillations 
are present, implying these modes are stable. As dis- 
cussed above, larger values for the vector and scalar field 
coupling parameters act to increase the size of the neg- 
ative region in the potential, which implies these cases 
are more likely to be unstable. We have extended our 
numerical analysis for extreme values of the parameter 
space, including variations of <p c , and find no evidence of 
any unstable modes. 

In figure 0] we plot the temporal evolution for the £ = 
1, 2 and 3 modes for K = k = 0.1 and ip c = 0.003. 
As expected there is no evidence of instabilities present 
for these, or any other values of the parameter space. 
We have further calculated evolutions for higher £'s and 
again find no evidence of any instabilities for this type of 
perturbation. 

Overlaid on figure 0] are damped sinusoidal oscillations 
where the frequency and damping times are those found 
using the WKB method in the following section. The 
WKB method is sensitive to the peak of the potential, 
whereas the time evolution accounts for the entire po- 
tential. One can see from these figures that the WKB 
method gives extremely accurate results, implying the 
negative region of the potential contributes to the evolu- 
tion of the perturbations on a level commensurate with 
the overall accuracy of the various schemes. 

We have further verified the stability of the solutions 
using the time independent form of the wave equation 
([331) . The method we used is based on the fact that, if 
unstable modes exist (wi m < 0), the boundary conditions 
for the radial perturbation function 'l'(r) become zero on 
both boundaries and the corresponding boundary value 
problem defined by equation ([35]) is self-adjoint pEi. lo6l]. 
This implies that, for the unstable modes, uj 2 is real and 
negative and to 2 should be greater than the minimum of 
the potential, V m - ln: for all of the bound states. Thus, 
in order to find unstable modes, we integrate equation 
([35]) with initial condition Xt'lr^-oo = for test values 
of lj 2 , starting from (J 2 = V m i n to w 2 = 0. An cigenfrc- 
quency is then found when the right boundary condition, 
^Ir.-xx) = 0, is fulfilled |67|- When applying this method 
to our problem, it turns out that no unstable modes ex- 



I4>l 
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FIG. 3: Typical temporal evolution of scalar field for £ — 
with K, = k = 0.1 and ip c — 0.003. The presence of oscilla- 
tions in the evolution implies non-zero oiR e which, due to the 
spherically symmetric nature of the background spacetime, 
implies the mode is stable. Due to the strong damping of 
the £ — mode, only a few oscillations can be seen in the 
evolution. 
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FIG. 4: Temporal evolution of scalar field for I — 1 (top 
panel), £ — 2 (middle) and £ = 3 (bottom) with K, = k — 0.1 
and ific = 0.003. The numerical time evolution is shown in 
black which includes the region of quasinormal ringing as well 
as the late-time tail. Overlaid in red is the damped sinu- 
soidal oscillations of the fundamental mode calculated using 
the WKB method (see section IVBl . 



ist for all of the studied values of the parameters and the 
black holes are therefore stable against the considered 
perturbations. 

B. Frequency domain calculations 

Given that no unstable modes exist, the next step is to 
calculate the QNM modes governed by equation (j35"[) . We 
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FIG. 5: The real (top) and the imaginary (bottom) part of 
the fundamental 1 = 2 QNM frequencies as a function of 
the scalar field coupling parameter k for several values of the 
vector field coupling parameter K. :— K + K+ — K4. The 
cosmological value of the scalar field is tp c = 0.003 and WR e 
and win, are shown in units of M. 
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FIG. 6: The real (top) and the imaginary (bottom) part of 
the fundamental 1 = 2 QNM frequencies as a function of the 
cosmological value of the scalar field tp c for several values of 
the parameters K. and k. Here, uiR e and oji m are shown in 
units of M 



have used two methods - a direct integration (shooting) 
method and also the WKB method. 

The first method we used is the shootin g m ethod in- 
troduced by Chandrasekhar and Detweiler [68[ (see also 
[69j and [66j|). A strong aspect of this method is that it 
takes into account the presence of the negative minimum 
of the potential, however it is also prone to suffer from 
numerical instabilities due to this minimum and also the 
fact that the first derivative of the potential ([55)1 is di- 
vergent at the left boundary, r = M/2. For this reason 
it was important to use an alternative method to con- 
firm our results. As can be seen in figure [T] and the top 
panel of figure [2J the negative minimum of the potential 
is extremely small compared to the positive maximum 
inrplyin g it is reasonable also to apply the WKB method 

In figures [5] and [6] we show the fundamental [n = 0) 
I = 2 QNM frequencies and damping times obtained from 
the WKB method for various values of the theory's cou- 
pling parameters 3 . As can be seen for relevant values of 
the parameters the frequencies only vary by a few percent 
throughout the range of the physically relevant param- 
eter space. We have displayed this in figures [5] and [5] 
only for the I — 2 mode however, as the associated er- 



I 




oj s /V2 




1 

2 
3 


0.0797 + i0.0749 
0.2071 + i0.0691 
0.3422 + i0.0684 
0.4779 + i0.0682 


0.0781 + i0.0742 
0.2071 + i0.0691 
0.3420 + i0.0684 
0.4776 + i0.0682 



3 For physically relevant values of the theory's parameters, the 
results obtained from the WKB and the shooting method differ 
by up to 1% and the difference becomes smaller for smaller values 
of the coupling parameters. 



TABLE I: Fundamental quasinormal modes for the first three 
spherical harmonics in the limit as the theory's various cou- 
pling parameters tend to zero, u) and for a Schwarzschild black 
hole in GR, u s - The Schwarzschild QNMs in GR are calcu- 
lated using the continued- fraction method, while u is calcu- 
lated using the shooting method which has an error in this 
limit of a few percent for the t = mode and significantly 
smaller for the higher modes. These numerical results con- 
firm the relation in equation (1291) . 



rors for the £ = and I = 1 modes are comparable to 
the change in frequencies associated with the variation of 
the parameter space. Nevertheless, our calculations have 
shown that the qualitative behaviour of the I = and 
1=1 modes are similar to the shown I = 2 case. 

Finally, from our calculations we have identified the 
limiting behaviour as the theory's various parameters go 
to zero which is shown in table HI The final column of ta- 
ble U shows the corresponding QNMs for a Schwarzschild 
black hole in GR. One can see that the errors produced 
are extremely small, even for the difficult to calculate 
I = case. These numerical results confirm the discus- 
sion of section llV Al and particularly the relation between 
GR and the limiting behaviour of TeVeS elucidated in 
equation ([2"fJ]) . 
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VI. CONCLUSIONS 



izcd to 



We have calculated here the stability and QNM spec- 
trum for spherically symmetric black holes in TcVcS and 
its generalization under perturbations of the scalar held 
intrinsic to the theory. Despite the existence of a small 
negative region in the potential of the wave equation, we 
have shown that these black holes are generically stable 
to such perturbations. In the limit as the various cou- 
pling parameters of the theory tend to zero, when one 
naively expects to recover GR, we have shown that the 
QNM spectrum tends to a value different to the scalar 
perturbations of Schwarzschild black holes in GR and in 
Brans-Dicke theory by a factor of \/2. The reason for this 
is clear: as one takes the limit of small coupling parame- 
ters the background vector field does not vanish. It is not 
until the parameters are set to zero that the background 
field vanishes, at which point the perturbation equation 
that we have explored is no longer a valid equation. 
In other words, this type of specific perturbation does 
not exist in General Relativity. 

We have explored herein only a small branch of exist- 
ing solutions in the strong-field regime in TeVeS. Indeed, 
even within the subset of spherically symmetric solutions 
with temporal pointing vector fields we have all but ig- 
nored a broad range of exotic solutions that involve either 
superluminal propagation of scalar perturbations, naked 
singularities or both. For these cases we have shown that 
the potential in the wave equation governing the pertur- 
bations generally (although there are some exceptions - 
see the Appendix) diverges at the critical horizon repre- 
senting either the naked singularity or the photon event 
horizon. The case that allows for superluminal pertur- 
bations is severely complicated by the fact that pertur- 
bations appearing within the photon horizon can escape 
to future null infinity. In fact, in these cases it is still 
not clear whether the central essential singularity (i.e. at 
r = 0) is naked to scalar perturbations in the sense that 
these perturbations at r = can evolve to future null 
infinity. We plan to explore the stability of these objects 
as well as those with naked singularities at the critical 
radius (i.e. those spacetimcs with a ^ 2) in future work. 



V 



2r 2 



M 
l + Yr 



-2(l+x) 



i-M 

2r 



-2(l-x) 



1(1+1) 



2r 



2r J 



2r 



2r, 
(A.l) 



where \ := a — 26±. Moreover, the coordinate, is now 
given by 



dr 



i+x 



dr \ 2r J \ 2r 



AI 



i-x 



(A.2) 



As with the more specific black hole case discussed in 
section HV CI the potential vanishes at spacelike infinity, 
however has the following limiting behaviour as one ap- 
proaches the horizon 



lim V(r) = 



for x > 2 
— oo for x < 2 



(A.3) 



It is known from section UTTI that 6 + > for all values of 
a, K. and k, while 6- < for certain values of a and K, 
depending on the relation given in (|19l) . 

Again, independent of the behaviour at r = M/2, the 
potential necessarily goes negative for some region of the 
spacetime, apart from the exceptional case of \ = 2. For 
the I = case, the potential has the following behaviour 
dependent on \: 



V(r; x < 2) < V 
V(r; x = 2) > V 
V(r; x > 2) < V 



M M 
ir<r<—, 
2 X 



M 



< r < 



X M 



(A.4) 
(A.5) 
(A.6) 



That is, for all cases excluding that where x = 2, the po- 
tential is negative for the immediate region of the space- 
time beyond r = M /2. At some finite value (either M/x 
for x < 2 or for x > 2) the potential is exactly 

zero, and is positive for all larger values of the radius. 
It is trivial to show that for all values of £ there exists 
a negative portion of the potential, however the mode 
that is most likely unstable is the fundamental I = 
mode. We leave the stability study of this broader class 
of spacctimes to future work. 



Appendix: The wave equation for naked singularities 
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In this appendix we show the more general wave equa- 
tion for arbitrary values of the metric parameter a, which 
is therefore relevant for perturbations of naked singular- 
ities amongst other things. Following the procedure out- 
lined in section IIV Al with an arbitrary a, leads to the 
wave equation pip , where the potential is now general- 
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